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ABSTRACT 

We wish to investigate the effects of cooling of the Compton cloud on the 
outflow formation rate in an accretion disk around a black hole. We carry out 
a time dependent numerical simulation where both the hydrodynamics and the 
radiative transfer processes are coupled together. We consider a two-component 
accretion flow in which the Keplerian disk is immersed into an accreting low- 
angular momentum flow (halo) around a black hole. The soft photons which 
originate from the Keplerian disk are inverse-Comptonized by the electrons in the 
halo and the region between the centrifugal pressure supported shocks and the 
horizon. We run several cases by changing the rate of the Keplerian disk and see 
the effects on the shock location and properties of the outflow and the spectrum. 
We show that as a result of Comptonization of the Compton cloud, the cloud 
becomes cooler with the increase in the Keplerian disk rate. As the resultant 
thermal pressure is reduced, the post-shock region collapses and the outflow rate 
is also reduced. Since the hard radiation is produced from the post-shock region, 
and the spectral slope increases with the reduction of the electron temperature, 
the cooling produces softer spectrum. We thus find a direct correlation between 
the spectral states and the outflow rates of an accreting black hole. 
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Introduction 



It is generally believed that the outflows and jets in a compact binary system containing 
black holes originate from the disk itself. There are several hydrodynamic al models of the for- 



matio n of outflows from the disks ra nging from the twin-exhaus t mo del of Bland ford fc Rees 



to the self-similar models of iBlandford fc Payne and IBlandford fc Begelman 



1999). Assuming that the outflows are transonic in nature, Fukue (jl983l ) and Chakrabarti 
( 119861 ) computed the velocity distribution without and with rotational motion in the flow 



and showed that the flow could become supersonic close to the black hole. Camenzind and 
his group extensively worked on the magnetize d jets and showed that th e acceleration and 



collimation of the jets could be achie yed (e.g.. lAppl fc Camenzind 



two component transonic flow model, IChakrabarti fc Titarchuk 1 fll995h pointed out that the 



19931 ). In a subsequent 



jets could be formed only from the inner part of the disk, which is also known as the Cen- 
trifugal pressure dominated boundary layer or CENBOL. The CENBOL is the Compton 
cloud in a two component flow and is the region between the accretion shock formed by 
centrifugal barrier at a few tens of Schwarzschild radii (r g ) and the inner sonic point located 
at ~ 2.5r 9 . If this region remains hot, it would emit hard radiation and the spectrum of the 
disk would be hard. The reverse is true if the Compton cloud (CENBOL) is cooled down 
and the spectrum would be soft. 

While the general picture of the outflow formation is thus understood and even corrob 



orate d by the radio observation of the base of the powerful jet, such as in M87 ( IJunor et al. 



19991 ) that the base of the jet is only a few tens of r g , a major question still remained: what 
fraction of the matter is driven out from the disk and what are the flow parameters on which 
this fra ction depends? In a numerical simulation using smoothed particle hydrodynamics 
(SPH), iMolteni et al. I (Il994[ ). showed that the outflow rates from an inviscid accretion flow 
strongly de pends on the outward centr ifug al force and 15-20 percent matt er can be driven out 
of the disk. IChakrabarti I ( 119981 . Il999h and lChattopadhyay et al. I (120041 ) estimated the ratio 
of the outflow rate to inflow rate analytically and found that the shock strength determines 
the ratio. For very strong and very weak shocks, the outflow rates are very small, while for 
the shocks of intermediate strength, the outflow rate is significant. These works were further 
refined by ISingh fc Chakrabarti I (120111 ) who studied the outflow rates self-consistently by 
modifying the Rankine-Hugoniot relations in presence of both energy and mass dissipation. 
These authors showed that with the increase of dissipation of ener gy in the flow, the ou tflow 
rate is greatly reduced. This is in line with the recent observations Fender et al. I ( 120101 ) that 
the spectrally soft states have less outflows. 



In the present paper, we concentrate on the numerical simulations of accretion flows 
around black holes which are coupled to radiative transfer. While computing the time vari- 
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ation of the velocity components, density and temperature we also compute the temporal 
dependence of the spectral properties. As a result, not only we compute the outflow proper- 
ties, we correlate them with the spectral properties - a first in this subject. Not surprisingly, 
we find that whenever the Compton cloud or the CENBOL is cooled down and the spectrum 
becomes softer, the flow, originating from CENBOL, losses its drive and the outflow rate is 
greatly reduced. 



Carrying out hydrodynamic simulations around black holes is not new. lHawley et al. 

(I1984T) pioneered the stu dy of accretion flows around black holes. Subsequent work of 



Hawley &: Balbus I ( 119921 ) with magnetized flow indicated how the angular momentum may 
be transported and m a tter m ay accrete efficiently These works used finite difference method. 



Igumenshchev et al. I ( 119961 ) studied two-dimensional flows but concentrated only the inner 



region of the disk, namely, the region less than 20r„. Their ma i n inte rest was to study 



the transonic nature of the flo w outside th e hori zon. iRyu et al. I (119951 ) studied unstable 



shocks in an inviscid flow, while IRyu et al. I ( 119971 ) studied the oscillations of the centrifugal 
pressure driven shocks when Rankine-Hugoniot relations are not satisfied. In a significant 
developm ent which connects th e time dependence of accretion flows with quasi-periodic os- 
cillations, iMolteni et al. I ( 119961 ) suggested that the shocks can oscillate whenever the cooling 
time scale roughly matches with the infall timescale in the post-shock region. Although only 
the bremsstrahlung cooling was used, this work provides a natural explan ation of the quasi- 



period ic oscillations or QPOs which are observed in black hole candidates. lLanzafame et al 



(119981 ) carried out the SPH simulation in a region with a radial extent of 50r g in two di- 
mensions and concentrated on the shock formation. They de monstrated that a high vis - 
cosity can remove th e shoc k from an accretion flow (see also, iGiri fc Chakrabarti 1120121 ). 
Igumenshchev et al. I (119981 ) used the finite difference method and allowed the heat gener- 
ated by viscosity to be radiated away or absorbed totally. The computational box was up to 
300r 9 , but the outer boundary condition was that of a near Keplerian flow having no radial 
velocity. The inner boundary was kept at 3r g . Thus the possibility of having a shock or 
the inner sonic point was excluded. The detailed radiative transfer, especially Comptoniza- 
tion was not included. To our knowledge, ours is a first step to include Comptonization 
coupled to hydrodynamic simulations where both the sources of soft photons (Keplerian 
disk) and the hot electrons (low angular momentum flows) are included. Evidence of such 
a two-component flow is present in many o f the observation s of the black hole candidates 



(ISmith et al. ll200lL 12002 ; ISoria et al. Il200ll ; IWu et al. 112002? ) . 



In the next Section, we discuss the geometry of the soft photon source and the Compton 
cloud in our Monte Carlo simulations. The variation of the thermodynamic quantities and 
other vital parameters are obtained inside the Keplerian disk and the Compton cloud which 
are required for the Monte Carlo simulations. In §3, we describe the simulation procedure 
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and in §4, we present the results of our simulations. Finally, in §5, we make concluding 
remarks. 



2. Geometry of the electron cloud and the soft photon source 

In Figure 1 , we present the cartoon diagram of our simulation set up for the Compton 
cloud with a specific angular momentum (angular momentum per unit mass) A = 1.73. The 
sub-Keplerian matter is injected from the outer boundary at R in = 100r 9 (r g = 2GM bh /c 2 ). 
The Keplerian disk resides at the equatorial plane of the cloud. The outer edge of this disk 
is located at R ou t = 200r 9 and it extends up to the marginally stable orbit R ms = 3r r 
At the centre, a black hole of mass M^u is located. The soft photons emerging out of the 
Keplerian disk are intercepted and reprocessed via Compton or inverse- Compton scattering 
by the sub-Keplerian matter. An injected photon may undergo no scattering at all or a 
single or multiple with the hot electrons in between its emergence from the Keplerian disk 
and its escape from the halo. The photons which enter the black hole are absorbed. 



2.1. Distribution of temperature and density inside the Compton cloud 



A realistic accretion disk is expected to be three-dimensional. Assuming axisymmetry. 
we have calculated the flow dynamics using a finite difference method which uses t he principle 
of total variation dimin ishing (TVD) to carry out hydrodynamic simulations (see. lRyu et al 



19971 ; iGiri et al. I l2010l . and references therein;). At each time step, we carry out Monte 
Carlo simulation to obtain the cooling/heating due to Comptonization. We incorporate the 
cooling/heating of each grid while executing the next time step of hydrodynamic simulation. 
The numerical simulation for the two-dimensional flow has been carried out with 512 x 512 
cells in a 100r g x 100r 9 box. We choose the units in a way that the outer boundary (Ri n ) 
is unity and the matter density at the outer boundary is also normalized to unity. We 
assume the black hole to b e non- rotating and we use the pseudo-Newtonian potential — 2 (r-i) 
( iPaczyhski fc Wiita 1119801 ) to calculate the flow geometry around a black hole (Here, r is in 



the unit of Schwarzschild radius r g .). Velocities and specific angular momenta are measured 
in units of c, the velocity of light and r g c respectively. In Figures (2) and (3) we show 
the snapshots of the density and temperature (in keV) profiles obtained in a steady state 
purely from the hydrodynamic simulation. For both the angular momenta cases, we find the 
presence of two shocks in the accretion flow in the steady sta te: one paraboloidal s hock and 
the other oblate spheroidal shock near the equatorial plane (IMolteni et al. Ill996l ). Due to 
the increased centrifugal barrier, the shock forms at a larger radius for the case A = 1.76. 
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2.2. Properties of the Keplerian disk 



In our simulation, the soft photons are produced from a Keplerian disk, the inner and 
the outer edges of which have been kept fixed at the marginally stable orbit R ms , and at 
Rout = 200r g respectively. The so urce of the soft photons h as a multi-color blackbody 



spectrum coming from a standard (jShakura fc Sunyaev I Il973l . hereafter SS73) Keplerian 
disk. The disk is assumed to be optically thick and the opacity due to free-free absorption 
is assumed to be more important than the opacity due to scattering. The emission is black 
body type with the local surface temperature (SS73): 

1/4 



T(r) « 5 x 10 7 (M 6/t )- 1 / 2 (M dl7 ) 1 / 4 (2r)- 3 / 4 



1 - 



K. 



The total number of photons emitted from the disk surface is obtained by integrating over 
all frequencies (u) and is given by, 



16tt 



fa 
he 



x 1.202057 



(T(r)f 



(2) 



The disk between radius r to r + 5r produces dN(r) number of soft photons. 

dN(r) = 47rr<5rif(r)77, 7 (r), (3) 
where, H(r) is the half height of the disk given by: 

/3 



H(r) = 10 5 M dl7 



1 - 



cm. 



(4) 



In the Monte Carlo simulation, we incorporated the directional effects of photons coming 
out of the Keplerian disk with the maximum number of photons emitted in the Z-direction 
and minimum number of photons are generated along the plane of the disk. Thus, in the 
absence of photon bending effects, the disk is invisible as seen edge on. The position of 
each emerging photon is randomized using the distribution function (Eq. 3). In the above 
equations, the mass of the black hole Mbh is measured in units of the mass of the Sun (M ), 
the disk accretion rate M^ 17 is in units of 10 1 7 gm/s. W e chose M^h = 10 in the rest of 
the paper. Generally, we follow iGhosh et al. I (120091 . 120101 ) while modeling the soft photon 
source. 



3. Simulation Procedure 



For a particular simulation, we use the Keplerian disk rate (rhd) and the sub-Keplerian 
halo rate (ihh) as parameters. The specific energy (e) and the specific angular momentum 
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(A) determines the hydrodynamics (shock location, number density and velocity variations 
etc.) and the thermal properties of the sub-Keplerian matter. We assume the absorbing 
boundary condition at r = 1.5 since any inward pointing photon at that radius would be 
sucked into the black hole. 



3.1. Details of the hydrodynamic simulation code 

To model the initial injection of matter, we consider an axisymmetric flow of gas in the 
pseudo-Newtonian gravitational field of a black hole of mass located at the centre in 
the cylindrical coordinates [R, 9, z]. We assume that at infinity, the gas pressure is negligible 
and the energy per unit mass vanishes. We also assume that the gravitational field of the 
black hole can be described by iPaczyhski fc Wiita I (119801 ) potential, 



GM, 



bh 



r — r n 



where, r = \/R 2 + z 2 . We assume a polytropic equation of state for the accreting (or, 
outflowing) matter, P = Kp" 1 , where, P and p are the isotropic pressure and the matter 
density respectively, 7 is the adiabatic index (assumed to be constant throughout the flow, 
and is related to the polytropic index n by 7 = 1 + 1/n) and K is related to the specific 
entropy of the flow s. K is not const ant but is allow e d to v ary du e to radiative proc esses. 
The details of the code is described in lMolteni et al. I (119961 ) and in lGiri et al. I (120101 ) . 



Our computational box occupies one quadrant of the R-z plane with < R < 100 and 
< z < 100. The incoming gas enters the box through the outer boundary, located at 
Ri n = 100. We have chosen the density of the incoming gas pi n = 1 for convenience. In 
the absence of self-gravity and cooling, the density is scaled out, rendering the simulation 
results valid for any accretion rate. As we are considering only energy flows while keeping 
the boundary of the numerical grid at a finite distance, we need the sound speed a (i.e., 
temperature) of the flow and the incoming velocity at the boundary points. In order to 
mimic the horizon of the black hole at the Schwarzschild radius, we place an absorbing inner 
boundary at r = 1.5r s , inside which all material is completely absorbed into the black hole. 
For the background matter (required to avoid division by zero) we use a stationary gas with 
density pb g = 10 -6 and sound speed (or temperature) the same as that of the incoming gas. 
Hence the incoming matter has a pressure 10 6 times larger than that of the background 
matter. All the calculations were performed with 512 x 512 cells, so each grid has a size of 
0.19 in units of the Schwarzschild radius. 
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3.2. Details of the radiative transfer code 



To begin a Monte-Carlo simulation, we generate photons from the Keplerian disk with 
randomized locations as mentioned in the earlier section. The energy of the soft photons at 
radiation temperature T(r) is calculated using the Planck's distribution formula, where the 



numb er density of photons (n^(E)) having an energy E is expressed by iPozdnyakov et al. 
jl983h . 



2C(3) 



b 3 E 2 (e bE 



1) 



(5) 



where, b = l/kT(r) and £(3) = Y1T ^ 3 = 1-202, the Riemann zeta function. Using another 
set of random numbers we obtained the direction of the injected photon and with yet another 
random number we obtained a target optical depth r c at which the scattering takes place. 
The photon was followed within the electron cloud till the optical depth (r) reached r c . The 
increase in optical depth (dr) during its traveling of a path of length dl inside the electron 
cloud is given by: dr = p n adl, where p n is the electron number density. 



The total scattering cross section a is given by Klein-Nishina formula: 



2vrr 2 



a 



x 



)ln{\ + x) + - + - 



x 



x 



2 (1 + x) 



where, x is given by, 



2E 



x 



7 1 - ji- 



(6) 



(7) 



e r/mc is the classical electron radius and m is the mass of the electron. 



We have assumed here that a photon of energy E and momentum ^ ft is scattered by an 

electron of energy / ymc 2 and momentum ~$ = 7m "v*", with 7 = ^1 — ^\ and p = O.v. 
At this point, a scattering is allowed to take place. The photon selects an electron and the 
energy exchange is computed using the Compton or inverse Compton scattering formula. 
The electrons are assumed to obey relativistic Maxwell distribution inside the Compton 
cloud. The number dN(p) of Maxwellian electrons having momentum between p to p + dp 
is expressed by, 

dN(p) oc exp[-(p 2 c 2 + m 2 c 4 ) 1/2 /kT e }dp. (8) 
In passing we wish to note that in the region of our interest (i.e. < 100r g ) and for stellar 



mass black holes, the free-free abso rption and/or emission is inefficient (e.g.. iNarayan fc Yi 
19951 ; iDas fc Chattopadhyay 1120081 ). This can be shown as follows: 



The Compton cooling rate for a t hermal distribution of nonre l ativistic electrons of num - 



ber density n e and temperature T e is ( IRybicki fc Lightman 111 9791 : iNiedzwiecki et al. 1119971 ): 
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C c = ^icn e aTU ra d, (U ra d is the radiation energy density, a? is the Thomson scattering 
cross-section, c is the speed of light, m e is the mass of the electron and k is the Boltzmann 
constant) while the bremsstrahlung cooling rate is given by C b = 1.4 x 10~ 27 n 2 e Te /2 . Using 
usual formulas for U ra d{r) and = 10, the ratio of the cooling rates turns out to be 

i? B = ^ = i.ixio 1 W- s (i-^)— ■ 

C b V r n e 

To compute R c , assume, rhh = 1.0 and m d = 1.0. For r ~ 100r 9 (pre-shock), T e ~ 10 s 
K, n e ~ 10 15 /cc, R c ~ 10 2 and for r ~ 10r g (post-shock), T e ~ 10 9 K, n e ~ 10 16 /cc 
R c ~ 10 4 . Thus both in the pre-shock and the post-shock regions, Comptonization is much 
more effective compared to bremsstrahlung. In the simulation we ignore bremsstrahlung. Of 
course, this conclusion depends on the mass of the black hole and should be considered for 
super-massive black holes, for instance. 



3.2.1. Calculation of energy reduction using Monte Carlo code: 

We divide the Keplerian disk in different annuli of width D(r) = 0.5. Each annulus 
having mean radius r is characterized by its average temperature T{r). The total number of 
photons emitted from the disk surface of each annulus can be calculated using Eq. 3. This 
total number comes out to be ~ ]^q41-42 p er seconc [ f or 7 j ld — \ q. In reality, one cannot 
inject these many number of photons in Monte Carlo simulation because of the limitation of 
computation time. So we replace this large number of photons by a lower number of bundles 
of photons, say, N comp (r) ~ 10 7 and calculate a weightage factor 

dN(r) 
fw ~N (r)' 

Clearly, from annulus to annulus, the number of photons in a bundle will vary. This is 
computed from the standard disk model and is used to compute the change of energy by 
Comptonization. When the injected photon is inverse-Comptonized (or, Comptonized) by an 
electron in a volume element of size dV, we assume that fw number of photons has suffered 
similar scattering with the electrons inside the volume element dV. If the energy loss (gain) 
per electron in this scattering is AE, we multiply this amount by fw and distribute this loss 
(gain) among all the electrons inside that particular volume element. This is continued for 
all the N comp (r) bundles of photons and the revised energy distribution is obtained. 
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3.2.2. Computation of the temperature distribution after cooling 

Since the hydrogen plasma considered here is ultra-relativistic (7 = | throughout the 
hydrodynamic simulation), thermal energy per particle is 3k bT where ks is Boltzmann 
constant, T is the temperature of the particle. The electrons are cooled by the inverse- 
Comptonization of the soft photons emitted from the Keplerian disk. The protons are cooled 
because of the Coulomb coupling with the electrons. Total number of electrons inside any 
box with the centre at location (ir, iz) is given by, 

dN e (ir,iz) = 4irrn e (ir,iz)drdz, (9) 

where, n e (ir,iz) is the electron number density at (ir,iz) location, and dr and dz represent 
the grid size along r and z directions respectively. So the total thermal energy in any box 
is given by 3kBT(ir,iz)dN e (ir,iz) = 127rrkBT(ir,iz)n e (ir,iz)drdz, where T(ir,iz) is the 
temperature at (ir, iz) grid. We calculate the total energy loss (gain) AE of electrons inside 
the box according to what is presented above and subtract that amount to get the new 
temperature of the electrons inside that box as 

A.E 

k B T new (ir,iz) = k B T i d (ir,iz) - — — . (10) 

3dN e {tr, iz) 



3.3. Coupling procedure 

The hydrodynami c and the radiative t ransfer codes are coupled together following the 



same procedure as in iGhosh et al. I (120111 ). Once a quasi steady state is achieved using 
the non-radiative hydro-code, we compute the radiation spectrum using the Monte Carlo 
code. This is the first approximation of the spectrum. To include the cooling in the coupled 
code, we follow these steps: (i) We calculate the velocity, density and temperature profiles 
of the electron cloud from the output of the hydro-code, (ii) Using the Monte Carlo code 
we calculate the spectrum, (iii) Electrons are cooled (heated up) by the inverse- Compton 
(Compton) scattering. We calculate the amount of heat loss (gain) by the electrons and its 
new temperature and energy distributions and (iv) Taking the new temperature and energy 
profiles as initial condition, we run the hydro-code for a period of time much shorter (by a 
factor of 20) than the cooling or infall time scale. Subsequently, we repeat the steps (i-iv). 
In this way, we see how the spectrum is modified and eventually converged as the iterations 
proceed. 

Earlier we mentioned that we injected N comp (r) bundles. Before the choice of N comp (r) ~ 10 
was made, we varied N comp (r) from ~ 10 3 to ~ 1.5 x 10 7 and conducted a series of runs 



7 



- 10 - 

to ensure that our final result does not suffer from any statistical effects. Figure 4 shows 
the average electron temperature of the cloud when N comp (r) was varied. There is a clear 
convergence in our result for N comp (r) > 3 x 10 6 . Thus, our choice of N comp (r) ~ 10 7 is 
quite safe. 

All the simulations are carried out assuming a stellar mass black hole (M bh = 1OM ). 
The procedures remain equally valid for massive/super-massive black holes though free-free 
absorption/emission may have to be included. We carry out the simulations for more than 
~ 10 dynamical time-scales. In reality, this corresponds to a few seconds in physical units 
for the chosen gridsize. 

4. Results and Discussions 



Table 1: Parameters used for the simulations. 



Case 


6, A 


m h 


m d 


la 


0.0021, 1.76 


1.0 


No Disk 


lb 


0.0021, 1.76 


1.0 


0.5 


Ic 


0.0021, 1.76 


1.0 


1.0 


Id 


0.0021, 1.76 


1.0 


2.0 


Ila 


0.0021, 1.73 


1.0 


No Disk 


lib 


0.0021, 1.73 


1.0 


0.5 


lie 


0.0021, 1.73 


1.0 


1.0 


lid 


0.0021, 1.73 


1.0 


2.0 



In Table 1, we list various Cases with all the simulation parameters used in the present 
paper. The specific energy (e) and specific angular momentum (A) of the sub-Keplerian halo 
are given in Column 2. Columns 3 and 4 give the halo ifn h ) and the disk ifn d ) accretion 
rates. The corresponding cases are marked in Column 1. In Cases la and Ila, no Keplerian 
disk was placed in the equatorial plane of the halo. These are non-radiative hydro dynamical 
simulations and no Compton cooing is included. To show the effects of Compton cooling on 
the hydrodynamics of the flow, the Cases I(b-d) and Il(b-d) are run for the same time as 
the Cases la and Ila. 



4.1. Properties of the shocks in presence of cooling 



In Figures 5(a) and 5(b), we present the time variation of the shock location (in units 
of r g ) for various Cases (marked on each curve) given in Table 1. All the solutions exhibit 
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oscillatory shocks. For no cooling, the higher angular momentum produces shocks at a higher 
radius, which is understandable, since the shock is primarily centrifugal force supported. 
However, as the cooling is increased the average shock location decreases since the cooling 
reduces the post-shock thermal pressure and the shock could not be sustained till higher 
thermal pressure is achieved at a smaller radius. The corresponding oscillations are also 
suppressed. The average shock location is found to be almost independent of the specific 
angular momentum at this stage. This is because the post-shock region is hotter in Cases 
I(a-d), and thus it cools at a very smaller time scale. In Figure 6, we show the colour map of 
the temperature distribution at the end of our simulation. We zoomed the region 50r 9 x 50r g . 
The specific angular momentum is 1.76 in the left panel and 1.73 in the right panel. Cases 
are marked. We note the collapse of the post-shock region as rhd is increased gradually. 
We take the post-shock region in each of these cases, and plot in Figures 7(a) and 7(b) the 
average temperatures of the post-shock region only for those cases where the cooling due to 
Comptonization is included. The average t emperature was obtained by the optical depth 
weighted averaging procedure prescribed in IChakrabarti fc Titarchuk I (119951 ). The average 
temperature in the post-shock region is reduced rapidly as the supply of the soft photons is 
increased. 



4.2. Effects of Comptonization on the outflow rates 

We now concentrate on how the outflow rate is affected by the Comptonization. Out- 
flows move to very large distances and thus must not be bound to the system, i.e., the 
specific energy should be positive. Matter should also be of higher entropy as it is likely 
to be relativistic. Because of this we wish to concentrate on the behaviour of matter which 
have highest energy and entropy. Though we injected matter at the outer edge with a con- 
stant specific energy, the energy of matter in the post-shock region is redistributed due to 
turbulence, Compton cooling and shock heating. Some entropy is generated as well. The 
high energy and high entropy matter escape in the form of a hollow cone around the axis. It 
is thus expected that if the post-shock region itself is collapsed due to Comptonization, the 
outflows will also be quenched. We show this effect in our result. In Figure 8, we present the 
specific energy distribution for both the specific angular momenta (left panel for A = 1.76 
and right panel for A = 1.73) for all the cooling Cases (marked in each box) at the end 
of our simulation. The velocity vectors are also plotted. The scale on the right gives the 
specific energy. First we note that the jets are stronger for higher angular momentum. This 
is because the post-shock region (between the shock and the inner sonic point close to the 
horizon) is hotter. Second, lesser and lesser amount of matter has higher energy as the 
cooling is increased. A similar observation could be made from Figure 9, where the entropy 
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distribution is plotted. The jet matter having upward pointing vectors have higher entropy. 
However, this region shrinks with the increase in Keplerian rate, as the cooling becomes 
significant the outward thermal drive is lost. Here, the velocity vectors are of length 0.05 at 
the outer boundary on the right and others are scaled accordingly. 

In order to quantify the decrease in the outflow rates with cooling, we define two types 
of outflow rates. One is M out which is defined to be the rate at which outward pointing flow 
leaves the computational grid. This will include both the high and low energy components 
of the flow. In Figures 10(a) and (b), we show the results of time variation of the ratio 
(Mout I Mi n ) for the four cases (marked in each box), Mj„ being the constant injection rate 
on the right boundary. While the ratio is clearly a time varying quantity, we observe that 
with the increase in cooling the ratio is dramatically reduced and indeed become almost 
saturated as soon as some cooling is introduced. Our rigorous finding once again verified 
what was long claimed to be the case, namely, the spectrally soft states (those having a 



Chakrabarti 


1999; 


Das et al. 


2001) 



Another measure of the outflow rate would be to concentrate only those matter which 
have high positive energy and high entropy. For concreteness, we concentrate only on matter 
outflowing within x = 20r g at the upper boundary of our computational grid. We define 
this to be Jm (= -vr^)- Figures 11(a) and 11(b) show time variation of J^. The different 
cases are marked on the curves. This outflow rate fluctuates with time. We easily find that 
the cooling process reduces this high energy component of matter drastically. Thus both 
the slow moving outflows and fast moving jets are affected by the Comptonization process 
at the base. 



4.3. Spectral properties of the Disk-Jet system 

In each simulation we also store the photons emerging out of the Computational grid 
after exchanging energy and momentum with the free electrons in the disk matter. When 
the Keplerian disk rate is increased, the injected soft photons go up, cooling every electron 
in the sub-Keplerian halo component. Thus, the relative availability of the soft photons and 
the hot electrons in the disk and the jet dictates whether the emergent photons would be 
spectrally soft or hard. In Figures 12(a) and 12(b), we show three spectra for each of the 
specific angular momenta: (a) A = 1.76 and (b) A = 1.73. The Cases are marked. We see 
that a spectrum is essentially made up of the soft bump (injected multi-colour black body 
spectrum from the Keplerian disk), and a Comptonized spectrum with an exponential cutoff 
- the cutoff energy being dictated by the electron cloud temperature. If we define the energy 
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spectral index a to be 1(E) oc E a in the region 1 — 10 keV, we note that a increases, i.e., 
the spectrum softens with the in crease in m&. This is consistent with the static model of 
Chakrabarti fc Titarchukl Jl995h . 



In reality, since the disk is not stationary, the spectrum also varies with time, and so is a. 
In Figure 13, we present the time variation of the spectral index for the different Cases. We 
draw the running mean through these variations. We clearly see that the spectral index goes 
up with the increase in the disk accretion rate. Thus, on an average, the spectrum softens. 
We also find that the spectrum oscillates quasi-periodically and the frequency is higher for 
higher cooling ra te. This agrees with the general observations that the QPO freq uency rises 
with luminosity flCui et al. Ill999l : ISobczak et ahlhoooi iMcClintock et al. 1120091 ). 



It is interesting to understand the physical reason behind the shock oscillations. I t has 



been argued by numerical simulation ( IMolteni et al. I Il996l ; IChakrabarti et al. 1 12004 ) this 



oscillation is due to a resonance effect. When the cooling time scale roughly agrees with the 
infall time scale, the post-shock region, i.e., the CENBOL may oscillate. As the CENBOL 
oscillates, the high energy part of the spectrum also oscillates since the CENBOL is mainly 
responsible to energize the soft photons coming out from the Keplerian disk to high energy. 
Indeed, when we compute the cooling time scale and the infall time scale, we find that in 
all the cases when oscillations occur, the ratio of cooling time scale and infall time scale on 
an average remains within forty percent of each other. Of course, in the absence of cooling 
effects, some ti mes the shock os cillations is possible when the Rankine-Hugoniot condition 
is not fulfilled (IRyu et al. Ill997l ). The amplitude of oscillation in this case is generally quite 
large and the result strongly depends on the specific angular momentum and specific energy 
of the flow. Observationally, different types of QPOs have been observed. This aspect 
requires a further study which is beyond the scope of this paper. 



5. Concluding remarks 

In this paper, we compute the effects of Compton cooling on the formation and rates of 
the outflows and jets. Most importantly, we study the variation of the spectral slope. We 
employ a hydrodynamic code based on the total variation diminishing method coupled to a 
Monte Carlo simulation to take care of the time dependent radiative transfer problem. We 
found many important results which can be briefly summarized as follows: (a) The increase 
in the Keplerian disk rate cools down the post-shock region of a sub-Keplerian flow and 
as a result the shock advances closer to the black hole, (b) The cooler and smaller sized 
post-shock region produced a weaker outflow/jet. (c) The spectrum becomes softer when 
the cooling is increased, (d) Quasi-periodic oscillations are produced as a result of near 
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resonance between the cooing time scale as the infall time scale. 



We coupled the hydrodynamics and Comptonization using a two step process. This 
is natural: Comptonization is strongly non-local process and the results are obtained by 
Monte Carlo method. However, the hydrodynamics code is a finite difference code. So at 
every instant both effects cannot be taken into account. The effects of radiative cooling or 
heating due to Comptonization is included in the hydrodynamics after each Monte Carlo 
step. After that the hydrodynamic code is advanced for a time period which is a factor of 
~ 20 lower than the cooling or hydrodynamic time scale, whichever is lesser before running 
the Monte-Carlo step once more. Thus our dynamics is affected by the radiative cooling or 
heating effects. 

In this context, it is important to note that the effects of radiation is not the same as 
in a single component flow. In a single component disk, such as the standard accretion disk, 
the radiation field is strong (especi ally when the accretion rate crosses the Eddington rate) 
as has been shown for example, by lOhsuga fc Mineshige I (120111 ). However, in our case, the 
radiation field is important only if it interacted with electrons. In a two component flow, 
there are two effects which affect the dynamics of the electrons (a) the thermal acceleration 
F t h depends on the gradient of thermal pressure and (b) the radiative acceleration F ra d 
depends on the number of Compton scattering and the amount of energy and momentum 
exchanged in the process (i.e., depends on local temperature). Both F t h and F ra d have strong 
dependence on both the accret ion rates. For a fixed halo ra te, F raA j, monotonically rises with 
disk rate but not exactly as in lOhsuga &: Mineshige I (120111 ). since the results depend on the 
optical depth of the halo. F t h need not rise monotonically, however. We find that with the 
rise of the disk rate, initial rising trend of F t h is reversed because the electrons are cooled 
down at a higher Keplerian disk rate, so the local temperature becomes low. We find that 
the ratio F ra d/F t h remains much less than unity even for ~ 10. 

Some of these results, as such, have been anticipated both theoretically and observa- 
tionally. From the theoretical point of view (IChakrabarti &: Titarchuk 1 119951). in the two 



component model showed that the spectrum becomes softer. Similarly, IChakrabarti I (Il999l ) 
found that the outflow rate becomes smaller when shock becomes weaker. The shock oscilla- 



tion model observations reported by iGallo et al. I (120041 ); iFender et al. I (120041 ) point to the 
same correlation between the spectral states and the outflow rates. However, it is the first 
time that we demonstrate the results with a complete coupled numerical simulations. The 
results we found are quite general and should be valid even for massive black holes, provided 
the cooling plays a role in the dynamics of the flow. 
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Fig. 1. — Schematic diagram of the geometry of our Monte Carlo simulations. Zigzag 
trajectories are the typical paths followed by the photons. The velocity vectors of the infalling 
matter inside the cloud are shown. The velocity vectors are plotted for A = 1.73. 




Fig. 2. — Logarithmic density (normalized unit) contours inside the halo for two different 
angular momenta (a) 1.76 and (b) 1.73. Logarithmic density contour levels are -4(0.25)2. 
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Fig. 3. — Temperature contours inside the halo. Velocity vectors are also provided. Here, 
temperatures are in keV. Temperature contours are 10(20)280. 




Fig. 4. — Variation of the average temperature (keV) of the post shock region with bundle of 
photons N comp for different Keplerian disk rates rhd, keeping the halo rate fixed at rhh = 1.0. 
Simulation cases (Table 1) are marked on each curve. Clearly, the temperature converges 
for N comp (r) > 3 x 10 6 . 
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Fig. 5. — The variation of shock location (in r g ) at the equatorial plane with time (in sec) 
for different Keplerian disk rates rhd, keeping the halo rate fixed at m h = 1.0. Simulation 
Cases are marked on each curve, (a) A = 1.76 and (b) A = 1.73. Cooling decreases the 
average shock location. 
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Fig. 6. — Colour map of final temperature distributions in the region (50r 9 x 50r g ) of the 
accretion disk for different disk rates are shown. The left panel is for A = 1.76 and the right 
is for A = 1.73. As is increased we find that, the high temperature region (dark zone, 
dark red online) shrinks. 




Fig. 7. — Variation of the average temperature (keV) of the post shock region with time (sec) 
for different Keplerian disk rates m^, keeping the halo rate fixed at = 1.0. Parameters 
are the same as in Figures 5a and 5b. 
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Fig. 8. — Colour map of final specific energy distribution inside the accretion disk for different 
disk rates. The high energy matter (dark zone, dark red online) are ejected outward as a 
hollow jet. The matter with a high energy flow decreases with the increase in disk rate. 
Velocity vectors at the injection boundary on the right is of length 0.05. 



Fig. 9. — Color map of the final entropy (K = -j^) distribution. Other parameters are as in 
Figure 8. The high entropy flow decreases as the disk rate increases. 
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Fig. 10. — Variations of 
and (b) A = 1.73. The Cases are marked in each panel. The outflow rate is the lowest for 
the highest Keplerian disk accretion rate (Cases are Id and lid). 
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Fig. 11. — Variations of J r - n (= 77^) with time for different is shown here. Here, M~ et 
and Mi n are the high entropy (also high energy) outflow and inflow rates, respectively. The 
left panel is for A = 1.76 and the right panel is for A = 1.73. The Cases are marked in each 
curve. 
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Fig. 12. — The final emitted spectra for different disk rates are shown for (a) A = 1.76 and 
(b) 1.73. Corresponding Cases are marked on each curve. The spectrum appears to become 
softer with the increase in m d . 
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Fig. 13. — Time variation of the spectral slope (a, 1(E) oc E~ a ) for different disk rates and 
their running averages (solid line, blue online) are shown. Different Cases are marked. We 
note that as the accretion rate goes up, the average a increases, i.e., the spectrum softens. 
The running mean gives an idea that there is a quasi-periodic variation of the spectral index 
- the frequency is higher for higher cooling rate. 



